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Chapter  10 


Permafrost  and  the 
Geothermal  Regimes 

A.  H.  Lachenbruch,  G.  W.  Greene, 
and  B.  V.  Marshall 


climatc  alT^V 13  ‘t/eUne 
active  climatic  change  that  has  been  Tprogretl 
throughout  the  past  century  has  increased  L 
annual  ground-surface  temperature  on  the  order  of 

permafrost  thickness  eIeSllyPwiU '^reduced  fl8nd 
1170  to  about  850  ft. 

ies  neat  the  shoreline  indicate  a  rapid  encroach 
men  of  the  Chukchi  Sea  several  thousand  years 
r  ‘'"f*  ,that  permafr08t  extends  under  the  marl 

at  !  H  fa  r -f  maJtlmum  distance  of  about  100  yards 
at  a  depth  of  200  or  300  ft.  y  8 

Preliminary  results  indicate  that  local  heat  flow 

ssr  “r°''  “*  ■—  - 


C,«  to  noL,..S"^.^*‘re 

sured  in  a  few  strateo-icaii,,  ,  ,  lasKa-  Temperatures  were  mea- 

selected  specimens  of  frozen  core  obtained^rom  ^  ^  ^ermal  proPerties 
Where  possible,  the  hole  locations ^  ™ 
important  role  played  by  bodies  of  wat  !»,,W  *  Permit  an  analysis  of  the 
simplified  by  the  uniformity  of  the  ^lthough  the  dermal  problems  are 

of  appreciable  heat ^  by  0)6  *■— 
information  on  thermal  conductivity  is  ne  Hod  u°f  °*  016  permafrost>  extensive 
terpretation  is  attempted  wHre  iust  hi  quantitative  in- 

ments  of  thermal  conductivity;  therefore  someth  °b*  h"  Iaboratory  measure- 
in  a  general  way  at  this  time.  6  results  can  be  reported  only 
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ST!" 01  “*  “r  «™»» 

surface  temperature  1.  also  beta  o“c  (aji ZrcWZ  i™7e£Sl 
temperatures  generally  obtain  to  depths  on  the  order  of  1000  ft  * 

ture  are  TeZed^f8  ^  ““  ° *  perermial!y  negative  tempera- 

that  miaht  hl  irrespective  of  the  physical  state  of  any  moisture 

nppnr  Present  (Muller,  1945).  Substantial  amounts  of  unfrozen  water  can 

SSI  bvTcfZ  r  *  r"UK  *P~«-  caused  4  Sec, Z 

capillarity  an^  ”  * 

f0rmPen^t°St  I"  b°Und,ed  b6l0W  by  a  °°C  isogeotherm,  which,  beneath  a  uni- 

wTc 

p  nds,  streams,  and  other  bodies  of  water  cause  large  local  thermal  anomalies 
which  propagate  downward  and  *  .  6  rindi  anomalies, 

of  permafrost  tZS,  temperature  and  hence  the  distribution 

-s^Sf  -  =  “-s  is 

AulTir'1  lrrepaarl,le8  “  ““  »•=  isogeothermal  ,urf.“e  “  ““ 

Jr-  -StfE* 

“  f  w/u,!h  >“•  luaxlmum  annual  temperature  1.  0"C.  pi,  practical 

srsviu™ be  mum  ™ •* 

into  "hrgeS.ln  °!f  temperature  at  the  ground  surface  propagate  slowly  downward 

r  rr 

-  “eS  ~  yl,ld  “» 


DETERMINATION  OF  NATURAL  EARTH  TEMPERATURES 

The  problem  of  determining  natural  earth  temperatures  has  two  parts 

oerZT^T  h°r  temperature  measurement  and  determination  of  the  tern- 
peragre,  disturbance  caused  by  drilling  ,h,  hole  in  which  temperature,  a”. 

Temperature  measurements  were  made  with  thermistors  spliced  into  multi- 
conductor  cables  at  regularly  spaced  intervals.  The  cables  were  lowered  into 
borehole,  and  left  in  place  for  at  least  1  year.  Periodic  me ZrmZ.  JZ 

brtdfe““L  ren'l  11’erl”1*t"  ”*d«  lrom  surface  with  a  Wheatstone 

bridge,  the  resistances  were  converted  to  temperatures  by  using  individual 

calibration  data  previously  determined  for  each  thermistor.  With  ^particular 
thermistors  and  measuring  equipment  used  in  this  study,  the  limit  of  accuracy  is 
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a  few  hundredths  of  a  centigrade  degree.  Where  possible,  cables  were  returned 
to  the  laboratory  for  recalibration  after  a  year  or  two  of  service  as  a  check  on 
the  stability  of  the  thermistors.  For  further  details  of  instrumentation  see 
Swartz  (1954),  Lachenbruch  and  Brewer  (1959),  and  Lachenbruch  et  al.  (1962). 

Thermal  cables  were  installed  in  four  holes,  designated  A,  B,  C,  and  D  (see 
Fig.  2),  to  depths  of  584,  1015,  1000,  and  1189  ft,  respectively.  Holes  A  and  B 
were  drilled  during  the  summer  of  1959  by  conventional  methods.  The  relatively 
warm  drilling  fluid  caused  thawing  and  sloughing  of  incompetent  wall  rock.  This 
thawing  and  sloughing  damaged  the  cables  and  rendered  many  of  the  thermistors, 
including  all  those  below  500  ft,  unreliable.  Holes  C  and  D  were  drilled  in  the 
summer  of  1960  by  the  Cold  Regions  Research  and  Engineering  Laboratory,  U.  S. 
Army  Corps  of  Engineers.  Their  drilling  with  refrigerated  fluid  (see  Chap.  6) 
required  special  techniques  developed  by  Robert  Lange,  who  directed  the  highly 
successful  drilling.  Drilling  with  refrigerated  fluid  had  three  great  advantages 
over  the  conventional  drilling  of  the  previous  year:  (1)  the  wall  was  maintained 
in  a  frozen  and  competent  condition,  and  thus  there  was  no  damage  to  the  cables; 
(2)  frozen  natural-state  permafrost  cores  were  obtained  for  laboratory  tests- 
and  (3)  relatively  little  thermal  disturbance  resulted;  so  natural  earth  tempera¬ 
tures  were  rapidly  restored. 

A  theoretical  analysis  of  the  dissipation  of  heat  introduced  by  drilling 
(Bullard,  1947;  Lachenbruch  and  Brewer,  1959)  shows  that  under  fairly  general 
conditions  the  return  of  the  temperature  O (a,/)  to  its  natural  predrilling  value 
©o(^),  at  any  depth  zt  can  be  described  by 

Q(i,i)  =  A  In  +  0OU)  ( l  »  s)  (1) 

where  t  is  the  time  since  the  drill  bit  reached  the  depth  z,  s  is  the  value  of  /  at 
the  cessation  of  drilling,  and  A  is  a  constant  sensitive  to  the  temperature  of  the 
drilling  fluid.  Thus,  a  semilogarithmic  plot  of  observed  temperature  vs.  (I  -  s)/t 
should  yield  a  straight  line  of  slope  A,  which,  when  extrapolated  to  (t  —  s)/l  =  1.0, 
i.e.,  infinite  time,  should  yield  the  natural,  or  equilibrium,  temperature  0^  The 
method  is  illustrated  in  Fig.  1  with  data  taken  near  the  300-ft  depth  in  each  of  the 
four  holes.  Note  the  rapid  approach  to  equilibrium,  indicated  by  the  small  slopes, 
of  the  two  holes,  C  and  D,  drilled  with  refrigerated  fluid.  The  positive  slope  for 
hole  C  indicates  that  the  permafrost  at  the  300-ft  depth  was  actually  cooled  by 
the  drilling  process.  The  equilibrium  temperature  is  determined  by  a  short 
extrapolation.  The  natural  earth  temperatures  determined  by  extrapolation  are 
presented  for  all  four  holes  in  Fig.  2  and  will  be  discussed  later. 


NEAR  SURFACE  TEMPERATURES 

The  air  temperature  near  the  ground  and  the  ground- surface  temperature  are 
subject  to  short-term  fluctuations  of  large  amplitude  corresponding  to  diurnal 
changes,  to  random  warm  spells  and  cold  snaps,  and  to  periodic  annual  seasonal 
changes.  Recording  three-pen  thermographs  installed  on  the  tundra  near  holes  C 
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Fi/t.  1 — Dissipation  uj  the  drilling  (list  n r  bailee  near  llic  .'iou-ll  depth.  Ilales  C 
and  I)  acre  drilled  nilli  re) rigeraled  fluid. 


and  D  permitted  a  study  of  the  downward  propagation  of  such  fluctuations.  These 
instruments  provided  a  continuous  record  of  the  temperature  (1)  in  the  air  (4  ft 
above  the  ground  surface  in  a  standard  U.  S.  Weather  Bureau  shelter),  (2)  1  in. 
beneath  the  ground  surface,  and  (3)  28  in.  beneath  the  ground  surface  near  the 
base  of  the  active  layer.  Thermograph  observations  near  hole  D  are  summarized 
on  the  upper  scale  in  Fig.  3.  Several  thermistors  installed  throughout  the  upper 
50  ft  of  holes  C  and  D  showed  the  short-term  temperature  fluctuations  in  the 
underlying  permafrost.  The  thermistors  were  read  weekly  or  biweekly  whenever 
possible.  Data  from  hole  D  are  presented  on  the  lower  temperature  scale  in 
Fig.  3. 

The  general  form  of  the  irregular  curves  obtained  from  air  and  from  the 
active  layer  was  obtained  by  smoothing  the  data  as  follows:  (1)  Mean  daily  tem¬ 
peratures  were  obtained  by  graphical  integration  of  the  thermograph  records; 
(2)  5-day  sliding  averages  were  computed  from  the  means;  and  (3)  these  were 
averaged  in  5-day  increments  and  plotted.  October  is  the  most  interesting  month 
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tig.  -J — Equilibrium  temperatures  in  holes  A,  B,  C.  anil  l).  S  hailing  indie  ales 
temperature  anomalies  associated  primarily  nilli  the  ocean,  the  lagoon,  or 
climatic  change. 


of  the  year  since  most  of  the  freeze-up  occurs  then.  Unfortunately  it  is  the  only 
month  for  which  active-layer  temperatures  are  not  presently  available.  Since 
more  near-surface  temperature  data  are  being  accumulated,  only  a  brief  qualita¬ 
tive  discussion  will  be  given  at  this  time. 

Perhaps  the  most  conspicuous  aspect  of  Fig.  3  is  the  smoothing  of  the  thermal 
fluctuations  as  they  penetrate  the  earth.  Even  when  they  are  artifically  smoothed 
and  plotted  on  a  contracted  temperature  scale,  the  surficial  temperatures  appear 
most  irregular  compared  with  the  unsmoothed  data  obtained  from  permafrost  a 
few  feet  below.  Most  of  this  smoothing  can  be  explained  by  the  fact  that  conducted 
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al  Me  D ",e  air  """ ac,il  c  la-™- M 


temperature  waves  are  attenuated  exponentially  with  depth,  the  negative  exponent 
being  proportional  to  the  square  root  of  the  frequency.  Thus  the  high-frequency 
components  (short-term  fluctuations)  of  the  surface  temperature  are  attenuated 
rapidly;  and  only  the  fundamental  frequency,  that  with  a  period  of  1  year  is 
significant  at  depths  greater  than  15  or  20  ft.  The  annual  wave  also  ultimately 
vanishes  as  its  amplitude  passes  from  about  15°C  at  the  surface  to  a  little  over 
C  at  - 15  ft  and  0.15  C  at  -50  ft.  At  depths  greater  than  70  ft,  the  annual  wave 
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is  no  longer  detectable  even  with  the  most  sensitive  thermal  elements.  System¬ 
atic  temperature  changes  that  persist  over  many  years,  climatic  changes,  for 
example  can,  of  course,  penetrate  to  much  greater  depths.  Such  a  change  is 
discussed  in  the  next  section. 

Figure  3  shows  that  temperature  fluctuations  at  the  ground  surface  are 
considerably  more  subdued  than  those  in  the  air  and  that  very  significant  attenua- 
tion  occurs  in  the  active  layer  between  the  surface  and  the  28-in.  depth.  Much 
of  this  attenuation  in  the  active  layer  is  the  result  of  seasonal  liberation  and 
absorption  of  latent  heat  from  the  freezing  and  thawing  of  interstitial  water  and 
is  not  related  to  the  normal  attenuation  in  conduction  discussed  above.  A  signifi¬ 
cant  shift  in  the  mean  annual  temperature  probably  occurs  between  the  air  and 
the  ground  surface,  but  this  cannot  be  investigated  until  more  data  are  available 
and  precise  calibration  checks  are  made  on  the  thermograph  elements 

Since  the  thermograph  element  at  -28  in.  is  within  a  few  inches  of  the  base 

°  vlre  1*yer>  11  can  be  marred  from  Fig.  3  that  the  freeze-up  of  the  active 

layer  at  this  site  was  completed  by  about  November  1  in  1960.  Very  rapid  cooling 
normally  follows  freeze-up  of  the  active  layer,  but  this  cooling  was  offset  some¬ 
what  at  the  28-in.  depth  by  the  effects  of  the  warm  spell  in  the  first  half  of 
November.  Rapid  cooling  is  conspicuous,  however,  at  the  5-  and  10-ft  depths  in 
eariy  December.  It  represents  a  rapid  adjustment  of  the  thermal  gradient  after 
the  retardation  of  winter  cooling  by  the  liberation  of  latent  heat  in  the  refreezing 
active  layer  (Lachenbruch  ei  at.,  1962),  6 

In  additl°"  t0  attenuation,  the  phase  retardation,  or  time  delay  in  the  down- 
r  propagation  of  temperature  fluctuations,  is  conspicuous  in  Fig.  3.  Both 
effects  are  well  illustrated  by  the  warm  spell  that  occurred  about  January  22  at 
:hSUr  --Ihe  6i'fects  of  the  warm  spell  appear  10  days  later  at  the  5-ft  depth 
f  u2,°  1  reductl0n  in  amplitude.  At  the  15-ft  depth  the  effects  are  scarcely 

detectable  as  a  slight  change  in  slope  in  mid-February.  Note  that,  because  of  the 
p  ase  retardation  in  the  annual  wave,  maximum  temperatures  occur  in  midwinter 
ana  minimum  temperatures  in  midsummer  at  the  30-ft  depth. 


RECENT  CLIMATIC  CHANGE 

temperatures  below  the  depths  of  annual  change  will  be  discussed  for  hole  D 
(Fig.  2)  first  because  this  hole  is  remote  from  bodies  of  water  and  their  compli¬ 
cating  thermal  effects.  The  almost  linear  portion  of  the  curve  below  -  500  ft 
represents  steady  heat  flow  from  the  earth  through  relatively  homogeneous 
ma  erials.  Extrapolation  of  this  straight  line  to  the  surface  (z  =  0)  yields  a 
temperature  of  -7.1°C.  The  results  of  this  extrapolation  strongly  suggest  that 
the  linear  temperature  profile  at  depth  developed  under  thermal  equilibrium  with 
a  mean  ground-surface  temperature  of  -7.1°C.  Extrapolation  of  the  measured 
profile  suggests  that  the  present  mean  surface  temperature  at  hole  D  is  about 
-5  C.  The  curvature  in  the  upper  500  ft  evidently  represents  a  departure  from 
thermal  equilibrium  caused  by  a  recent  trend  toward  higher  mean  ground-surface 
temperatures.  The  nature  of  this  trend,  which  presumably  represents  a  climatic 
change,  can  be  investigated  with  the  theory  of  heat  conduction.  The  problem  has 
been  discussed  by  Birch  (1948). 
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perafure  wet^calcu  tT  hm  ^  the  Changing  mean  annual  surface  tem- 
to  have  started  to  orSl  fi  ™agnitude  o£  the  change  and  the  time  it  would  have 
observed  proSe  af  hrfeD °"  P6rmafr°St  ^-lures.  If  the 
anomaly  can  be  determined  (cur™ D fJ“  toe  cUmrtic 

temperature,  e,,  was  assumed  to  have  varied  as  '  3nnUal  surface 


Os  = 


6. 


0  £  A  S  /* 


(2) 


where  A  represents  time  before  the  present  and  s  io 

T-1VC  The  d°f  ‘I6  The  ValUe  determined  fZ^r'veD 

is  best  represented  bv  !  u'  ,T  3  h3d  a  ma«nitude  of  2°C.  If  the  change 
started  about  ?nn  y  ,ncreasi«g  surface  temperature  //  1  it 

m  t:  ?  “ a  -r-  of  “  2,^‘ 

about  125  and  150  years  ago,  respectively.' The  four  models  J'T  i°  Tand^rt 
Z  IZIZTZZZ  o’o2°C  TheP  n‘  error’ o/’a^tO.orc 

mentr^Ts^Alor^le1*8'  ?  T  °f  u-ertainty  of  the  measure- 
teU  the  same  genera,  lr  -  the  ,  f  t"  "*  m°delS-  Nevertheless.  they  all 
last  hundred  years, f  and  it' increased  th  was  essentially  an  event  of  the 

hole  D  by  2“  to  2«/*C A  siZZ  h  ^  ^ound-surface  temperature  at 

-  — 

all^ge^e^  temperature  are  Conspicuous 

(Misener,  1955)  B3y’  Northwest  Territories,  Canada 


tes^we^ttypicarand  <j’a‘  lhc  early  samples 

7  x  to  *  cal/°C/sec/cm.  This  yields  i  heat  flow  of  ihn  \  ty /°’e  * "■’  f°rmatton  is  about 
coincidence,  is  still  “close  to*  he  worldwide  X  10  c;,1/cni2/sec.  which,  by 

interim  have  resulted  in  an  upward  revision  of*  the  ,SinCn  many  de,el'minations  in  the 
Is  7  instead  of  5,5  *  10~3  cal/cm/’C/sec  as  oriirtn.llv  g°'  ®eCaa8e  ,ht'  mo'ln  conductivity 
in  calculation  of  shoreline  movents  and  ciimaHc  ehan^  ’  diffush  ity  values  used 
tlu'  time  determined  from  these  analyses  should  h.  ‘  "8e  we| 0  Proportionally  low:  hence 
general  results  are  unchanged  *  6  *"  r*vtoed  dow"'™'d  by  25<'c  or  so.  The 
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Thus,  an  analysis  of  the  case  n  =  V,  (Fitr  4  K\ 

warming  is  equivalent  to  the  introdimtin  6 *  ,!  P  rt  b  ’  reveals  0131  016  climatic 
the  rate  of  Lut  20  cal/^m  ‘  ~  ,  ^at  through  the  ground  surface  at 

about  half  as  great  as  £  steadvynar  ^k^  PaSt  eight  decades-  rate  is 

•Hi.  only  ,  "oi'Sjl Z  H,7t  £??  T  eM'"‘ 

at  temperatures  above  the  mean  and  leaves  it  ^  the  ground  surface 

with  the  changing  seasons.  1  temperatures  below  the  mean 

r.cJ.^C PT'  !>el™  -500  * 

...  iie.! 

cannot  be  accounted  for  as  nart  nt  ,  /  ,  mat  ““  recent  warming  trend 

mean.  The  data  imply  further  that  llrg^swle^ysfeSi^ch  ***  h  l0ng'term 
curred  in  the  last  few  thousand  years.  ?  systematic  changes  have  not  oc- 

thouslfyearT^  TTn  ™  ‘°  P6rSiSt  for  *-eral 

configuration  represented  by  a  straight’ Une^^^ 
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depth  and  intersecting  the  surface  at  about  -5°C.  Under  such  conditions  the  0°C 
isotherm  (bottom  of  permafrost)  would  lie  at  a  depth  of  only  about  -850  ft.  In 
this  sense  about  25%  of  the  permafrost  beneath  Ogotoruk  Valley  is  the  product  of 
an  extinct  climate. 

It  is  important  to  emphasize  that  the  increase  in  mean  annual  ground-surface 
temperature  does  not  necessarily  imply  that  the  mean  annual  air  temperature  has 
increased.  Subtle  changes  in  other  climatic  parameters  could  account  for  the 
trend.  Although  this  change  in  mean  ground-surface  temperature  might  not,  in 
itself,  be  important  to  local  biological  systems,  it  probably  was  accompanied  by 
marked  changes  in  other  more  biologically  significant  parameters,  such  as 
amplitude  of  seasonal  temperature  variation  or  quantity  of  winter  snow  or  summer 
rain.  The  mean  annual  ground-surface  temperature  just  happens  to  be  the  quantity 
accessible  to  a  geothermal  analysis. 


HOLE  B  AND  THE  SURFACE  MICROENVIRONMENT 


The  thermal  profile  at  hole  B  has  a  similar  but  less  pronounced  curvature 
than  that  at  hole  D  (Fig.  2).  Since  only  three  uncertain  measurements  are  avail¬ 
able  below  the  400-ft  depth  in  hole  B,  it  is  not  possible  to  determine  tne  gradient 
accurately  at  depth.  One  reasonable  interpretation  is  represented  by  the  straight 
line  drawn  asymptotic  to  the  profile  for  hole  B  in  Fig.  2.  This  interpretation 
would  imply  that  before  the  recent  climatic  change  the  mean  surface  temperatures 
at  holes  B  and  D  were  the  same,  or -7.1°C.  The  corresponding  geothermal 
anomaly  is  plotted  as  curve  B  in  Fig.  4,  part  a.  Curve  B  can  be  obtained  within 
instrumental  error  by  multiplying  the  abscissa  for  curve  D  by  a  factor  of  0.62. 
Thus  the  five  models  represented  in  Fig.  4,  part  b,  account  for  curve  B  within 
instrumental  error  if  the  temperature  scale  is  reduced  by  0.62.  Hence  a  total 
change  of  mean  surface  temperature  at  holeDof  2°C  corresponds  to  one  at  hole  B 
of  about  1.25°C. 

The  difference  in  the  present  mean  annual  surface  temperatures  over  the  half 
mile  or  so  between  holes  D  and  B  seems  to  be  real  although  the  value  of  this 
difference,  about  V/C,  can  be  adjusted  somewhat,  depending  on  the  model  used  as 
a  basis  for  interpretation.  The  evidence  is  inconclusive,  but  it  does  suggest  that 
the  mean  annual  ground-surface  temperatures  at  the  two  sites  were  not  signifi¬ 
cantly  different  before  the  climatic  change.  Thus  it  is  likely  that  the  microen¬ 
vironmental  difference  causing  the  present  disparity  between  the  surface  tem¬ 
peratures  at  holes  B  and  D  probably  evolved  during  the  last  century  and  might 
well  be  in  the  process  of  change  today.  As  an  example,  the  surface  in  the  vicinity 
of  hole  D  has  high  microrelief  with  Eriofihorum  tussocks  interspersed  with  bare 
mineral  soil.  At  hole  B  the  microrelief  is  less  extreme,  and  the  vegetal  mat  is 
more  continuous.  If  the  rough  surface  near  hole  D  traps  more  drifting  snow,  this 
could  account  for  the  difference  (Lachenbruch,  1959).  The  existing  data  could  then 
be  interpreted  as  an  indication  that  the  extreme  microrelief  developed  recently 
near  hole  D,  perhaps  as  a  result  of  the  climatic  change.  Several  other  possible 
causes  could  be  mentioned,  but  the  identification  of  the  most  likely  factors  is  best 
left  to  the  ecologists. 
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(insert  LV’ it  VII  dnUed  at  Varying  distances  ^om  ‘he  shoreline 

between  ^  by  the  Shading  in  Fl^  2  that  the  major  difference 

and  u* L  o  ?!  !  Pr°£lles  is  the  result  of  thermal  effects  of  the  ocean 
at  hS  '.  r  bodies  Of  water  can  cause  surprisingly  large  geothermal  effects 

temneL  kS  b6CaUSe  the  Perennially  unfrozen  bottom  sediments  have  a  mean 
temperature  above  the  freezing  temperature  of  the  water  that  contacts  teem  bS 

s  xTJTo  sfce  ‘v'ner*uy  m“y  a'8"*s  « c 

earth  s  c  ,  °r  when  a  shoreline  moves,  the  temperature  of  a  portion  of  the 
gS  d0~7dCa  dTt  ranCaUy-  The  resolting  geothermal  disturbance  propa- 
r‘ ,hC  'Srt6;  *  *a^tes  thousands  of  years  before 

P  .  .  Unding  emer?ent  surface  (assumed  uniform) 

1961) 1 'Z  Z  T°m  °f  Z  ChUkChi  Sea  Sl0pes  very  «ently  (Scholl  and  Sainsbury 
5“  cXcT  ,™”,arPS,’  11,5  ‘““““J  “»  estimated  (ioS 
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pointBrert  haS  Sh°Wn  that  the  temPerature  of  shallow  coastal  waters  near 

S  al  0  eTo"o  ‘7J‘Z  “  “«»—“»«>  tanSrJS 
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communication,  1960)  sugge.ts  lhat  ISM  mean  water 
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lS  °'75  C<  The  corresponding  averages  for  Point  Barrow  and 


+M.  C.  Brewer, 
published. 


Ice.  Water,  and  Bottom  Temperatures  near  Point  Barrow.  Alaska,  u„- 


|t|l  i  i 


160 


LACHENBRUCH,  GREENE,  AND  MARSHALL 


JScJS  f,rt0m  ul  data  °f  Bre""r  and  B100m  are  ab0Ut and  1 0.90  C, 
respectively  It  will  be  assumed  that  the  mean  bottom-sediment  temperature  has 
remained  close  to  this  valup  f+n  r7R°r,\  <■,-»  4-k 

Thi«  *  (  0,75  C)  for  the  Past  several  thousand  years. 

This  assumption  is  supported  by  evidence  from  Point  Barrow  that  the  recent 

lS  and  bTevii^  '‘T  6ffeCt  °"  SUbmereed  surfaces  (Lachenbruch  dal., 
1962)  and  by  evidence  from  hole  D  that  indicates  a  uniform  gradient  at  depth 

Srew  l„SUy  “sSyS,ema,‘C  *«*»«» 


at  hole  C  (386  ft  H  J  >  fSUmptlons'  the  the™ial  anomaly  caused  by  the  sea 

bruch  f9  7!  in  f  f(  ee"  C°mpUteCl  fr0m  heat-conduction  theory  (Lachen¬ 
bruch,  1957a).  In  Fig.  5  the  curve  marked  /  c  represents  the  effect  the  sea 

would  have  if  the  shoreline  had  been  in  its 


temperature  anomaly,  °c 


present  position  long  enough  for  thermal  equi¬ 
librium  to  obtain.  The  other  two  broken  curves 
in  Fig,  5  represent  the  thermal  effect  to  be 
expected  at  hole  C  if  the  shoreline  had  moved 
suddenly  to  its  present  position  4000  or  1000 
years  ago. 

Curve  C-D  (Fig.  5)  represents  the  differ¬ 
ence  between  the  temperatures  measured  at 
holes  C  and  D.  It  can  be  shown  that  hole  D  is 
so  far  inland  (3900  ft)  that  temperatures  there 
are  unaffected  by  the  ocean.  Hence,  if  other 
conditions  are  equal  at  the  two  sites,  curve 
C-D  represents  tiie  thermal  disturbance  caused 
by  the  sea  at  hole  C.  Curve  C-B  represents  the 
difference  between  measured  temperatures  at 
holes  C  and  B  (about  2000  ft  inland).  From  the 
shapes  of  curves  C-B  and  C-D,  it  is  clear  that 
the  theoretical  model  (broken  lines)  is  too 
simple  to  account  for  details  of  the  ocean 
disturbance.  Since  C-D  does  not  increase  with 
depth,  however,  the  present  shoreline  position 
evidently  has  not  been  occupied  long  enough  for 
the  ground  to  attain  thermal  equilibrium  with 
the  sea.  The  simple  model  of  a  rapid  shoreline 
transgression  4000  years  ago  fits  curve  C-D 
very  well  below  700  ft  and  hence  agrees  with 
the  measured  depth  of  permafrost  at  hole  C 
(945  ft).  Near  the  surface.  C-B  and  C-D  bracket 


iiH.  ■>—. Analysis  „J  thermal  meets 
Lilli)'  C  -it  represents  tempi  rtth, re  iliiiereiiees  he- 
/l"'i :"‘!rs  c  ""<<  «•  Cure,  C-D  represents  tem¬ 
pera  fan  dntc  fences  he! m  en  holes  C  mi, I  D.  linden 
hues  represent  the  computed  cure/  at  hate  C  01  a 
rapid  shoreline  transMressimi  In  the  pres,  nt  position 
1  years  ayo. 
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the  theoretical  curves;  thus  their  departure  from  the  model  can  be  explained  if 
the  mean  surface  temperature  at  hole  C  is  intermediate  between  that  at  hole  B 
and  at  hole  D.  The  departure  of  C-B  from  C-D  at  depth  is  evidently  the  result  of 
local  topographic  effects  at  hole  B  and  possibly  of  a  minor  difference  in  thermal 
conductivity. 

If  the  departure  at  intermediate  depth  between  curve  C-D  and  the  curve 
!  4000  were  to  be  explained  in  terms  of  shoreline  history,  we  would  be  forced 

to  infer  that  several  thousand  years  ago  the  transgressing  shoreline  "overshot" 
by  several  hundred  feet  and  then  retreated  to  its  present  position  in  the  last 
century.  However,  such  an  explanation  is  precluded  by  the  present  beach  mor¬ 
phology.  Although  alternative  explanations  are  possible,  the  most  likely  explana¬ 
tion  is  that  the  disparity  represents  the  combined  effect  of  a  low- conductivity 
layer  in  hole  C  between-200  and  -300  ft,  as  indicated  by  the  change  in  gradient 
(Fig.  2),  and  the  neglected  effects  of  the  latent  heat  of  melting  interstitial  ice  at 
the  rising  base  of  permafrost.  A  preliminary  calculation  accounting  for  these 
effects  suggests  that  the  thermal  data  are  consistent  with  a  rapid  shoreline 
transgression  occurring  about  6000  years  ago.t  Although  this  date  is  consistent 
with  that  inferred  by  Giddings  (1960)  from  archeological  studies  and  by  Moore 
(1960)  from  study  of  the  beaches,  it  could  be  revised  upward  or  downward  by  as 
much  as  50'<  when  additional  thermal  conductivity  values  become  available  and 
make  more  refined  calculations  possible. 

The  distribution  of  permafrost  corresponding  to  the  theoretical  models 
represented  in  Fig.  5  is  shown  in  Fig.  6.  The  theoretical  curve  corresponding 
to  a  rapid  transgression  4000  years  ago  gives  the  depth  of  permafrost  correctly 
at  the  two  places  where  it  was  measured  (holes  C  and  D)  and  probably  provides 
a  reasonable  approximation  to  the  actual  permafrost  distribution  even  though  the 
actual  date  of  transgression  is  underestimated  because  of  the  neglected  effects 


t'i^.  PennaJ rosl  distribution  computed  for  a  shoreline  transgression  t 

years  ago.  The  c/Jects  oj  latent  heat  have  been  neglected.  Curve  l  =  loon  years 
is  close  to  present  conditions  although  the  time  since  transgression  is  probably 
greater.  Asterisks  denote  the  measured  depth  oj  permafrost. 


tSee  footnote  on  p.  15(i. 
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twiTodercJt  r°Utby/hiCh  timC  °f  movement  ^  underestimated  by 

between  the  curves  /  I  0ld  /  motJT  *°  meW  0,6  interstitial  ice 

to  be  about  1%  or  less.  '  aVerage  lce  content  is  believed 

As  implied  in  Fig.  2,  the  major  cause  of  the  difference  between  the  tem 

rcrr c  ls  “ua  « ,te  ■*«»■ « 

nne  trend  is  projected  across  the  mouth  of  the  lagoon,  holes  A  and  C  are  about 

bf  about6  th  06  ^°m  U;  h6nCe  the  thermal  effects  of  the  sea  a«  expectedto 
...  e  same  in  each  hole,  and  the  convergence  at  depth  of  the  temperature 

dP'01'  r  r-  A  “  0  <F‘S-  !>  “  «*  ‘»rpr.,lns.  Since  .he  I'ZZZ 
itl  tH  ,  »°°n  Vary  appreciably  {rom  y^r  to  year,  it  is  not  possible  to  study 
its  thermal  effect  in  detail  by  theoretical  methods.  Nevertheless,  application  of 
the  three-dimensional  heat-conduction  model  (Lachenbruch,  1957b)  shows  that  the 
difference  between  the  temperature  profiles  at  holes  A  a^  C  2  of  the  e^ecSd 

g  neral  form  and  magnitude  if  the  lagoon  is  assumed  to  have  occupied  its  present 
position  for  the  last  few  thousand  years.  occupied  its  present 

PRELIMINARY  ESTIMATE  OF  EARTH  HEAT  FLOW 

0.  JSZc ln  T D  *” ”• 

identified  lith  th  t  temperature  with  depth  below  -  500  ft  can  be 

identiRed  with  the  outward  flow  of  heat  from  the  earth’s  interior.  The  linearity 

of  the  profile  suggests  that  the  rocks  have  a  relatively  uniform  thermal  con¬ 
ductivity  when  viewed  on  the  scale  of  the  thermistor  spacing.  Although  detailed 
information  is  not  yet  available  on  the  conductivity,  we  ran  about24  tests  on 
frozen  core  selected  at  random  while  we  were  developing  and  Jefinii  the  con- 
^uctivlty  apfiaratus.  The  core  samples  yielded  an  average  value  0f5  5  x  10"3 

low7  500  ft  CC'  h,  nPlmt  thiS  iV6rage  Value  by  tbe  mean  thermal  gradient  be¬ 
low  -500  ft  in  hole  D  (19.8  x  lO^C/cm)  yields  a  heat  flow  of  about  1.1  x  HT6 

vrllr,/366’  ?iC!“  18  Ci°Se  t0  1116  worldwide  average  as  we  know  it  today.!  The 
greatest  uncertainty  in  this  calculation  lies  in  the  selection  of  the  mean  value  of 

morl  T  COnductmty-  which  miebt  be  adjusted  by  as  much  as  10  or  15%  after 
more  extensive  measurements  are  made.  The  neglected  thermal  effects  of 
topographic  relief  would  affect  the  results  by  only  about  2%,  and  the  effects  of  a 
ypothetical  change  in  mean  ground-surface  temperature  at  the  end  of  Wisconsin 
gbacution  by  possibly  10%.  Since  the  area  was  not  glaciated  in ^scoMin  timf 

SionJl’h  fY*  n°  evidence  that  any  such  change  occurred;  nor  is  there  any 

TZ  "?cool^PeCUlating  “  ^  a  hyP°thetiCal  Cha-e  a 

19„feat  ®tudies  at  Resolute  Bay.  Northwest  Territories,  Canada  (Misener 

s,g  “Srr;  and  LerTOX’  1962)  have  led  aonie  authors  to 
2  1,7  ,  1  America  is  a  region  of  greater  than  normal  heat  flow 

he  results  from  Ogotoruk  Valley  are  not  consistent  with  such  a  generalization 


tSee  footnote  on  p.  156. 
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